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Abstract 

The possible existence of light sterile neutrinos is of great interest in many areas of particle 
physics, astrophysics, and cosmology. Furthermore, should the MiniBooNE experiment at 
Fermilab confirm the LSND oscillation signal, then new measurements are required to identify 
the mechanism responsible for these oscillations. Possibilities include sterile neutrinos, CP 
or CPT violation, variable mass neutrinos, Lorentz violation, and extra dimensions. In this 
paper, we consider an experiment at a stopped pion neutrino source to determine if active- 
sterile neutrino oscillations with Am 2 greater than 0.1 eV 2 can account for the signal. By 
exploiting stopped tt + decay to produce a monoenergetic source, and measuring the rate of 
the neutral current reaction v x 12 C — ► v x 12 C*(15.11) as a function of distance from the source, 
we show that a convincing test for active-sterile neutrino oscillations can be performed. 



1 Introduction 

The understanding of neutrino mass has undergone a revolution over the last ten years. Neu- 
trino mass, via oscillations, is now experimentally established pQ. The Standard Model may 
now need to accommodate right handed neutrinos, vr, that are sterile with respect to the weak 
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interaction. The new issue to be addressed concerns the mass spectrum of these sterile neu- 
trinos. Although many theorists postulate that i/r are heavy via the see-saw mechanism (2], 
there is no experimental evidence for the mass spectrum of sterile neutrinos or for the number 
of generations 0. In fact, light sterile neutrinos have been invoked to explain many puzzles, 
such as the LSND result |H, the solar neutrino energy spectrum anomaly [3], the R-process 
in type II supernovae [B], pulsar kicks [7j, dark matter [S], the perplexing dark energy and 
most recently, extra dimensions |1U| . Thus, much experimental and theoretical justification 
exists to motivate a search for light sterile neutrinos above the mass of 0.1 eV 2 . 

The LSND result points to active-sterile neutrino oscillations in the Am 2 range of 0.1 
to 2eV 2 A recent analysis of a (3+2) sterile neutrino model indicates the possibility 

of a second sterile neutrino solution around 20 eV 2 0]. In addition, a supernovae R-process 
analysis favors a 3 < Am 2 < 70 eV 2 sterile neutrino [Hj. We propose an experiment that is 
capable of the direct detection of active-sterile neutrino oscillations in this mass region. If 
we interpret the LSND result as an intermediate sterile state, i.e. — > P B — ► P e , where the 
observed LSND oscillation probability is P(P^ — > v e ) = (0.264 ± 0.081)% then coupling 
this with short baseline reactor experimental limits on v e disappearance of -P(z^ e — > u B ) < 10% 
[T2] . we can predict the disappearance to sterile neutrinos to be Piy^ — > v s ) > 10%. The 
above arguments set the scale for the Am 2 and sin 2 29 sensitivity that needs to be achieved. 

2 Beyond the Standard Model 

Given the prospect of a new intense stopped pion neutrino source, together with a posi- 
tive signal from MiniBooNE a next generation short baseline neutrino experiment needs 
to be carefully considered. A neutrino detector at the Spallation Neutron Source (SNS) 
has already been discussed |14| that demonstrates how an improved measurement in the 
LSND/MiniBooNE region can be made. Improving the Am 2 and sin 2 29 measurement of 
any positive signal from MiniBooNE would be important. However, the nagging question 
would still remain as to the new physics behind the three distinct Am 2 regions that have been 
measured with solar, reactor, atmospheric, and short baseline detectors. Several possibilities 
at present exist, such as sterile neutrinos [I], CP or CPT violation J3]> m ass varying neutrinos 
[5], Lorentz violation [T2|, and extra dimensions jlOj . 

The CP or CPT violation possibility could be readily tested with a MiniBooNE/LSND 
style detector coupled with a short duty-factor stopped pion beam which allows, via timing 
cuts, a separation of Vn and Vn events. Thus, a direct comparison of oscillation probabilities 
from these two simultaneous data sets will allow a test of CPT and CP violation models. Mass 
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varying neutrinos can be tested by changing the density of material between the source and 
detector. The Lorentz violation model can be tested by looking for sidereal variations in the 
neutrino flux. These three tests will be discussed in more detail in a future paper. 

This paper discusses the first possibility, that of sterile neutrinos. To demonstrate active- 
sterile neutrino oscillations directly, and achieve the best sensitivity, requires a two detector 
setup and a stopped pion neutrino source. The following sections describe how this can be 
done. 

3 Intense Stopped Pion Neutrino Sources 

To detect sterile neutrino effects in the mass range of LSND and above requires an intense 
source of well characterized neutrinos. Stopped pion decay from a low energy proton beam is 
such a source. Neutrinos from stopped pion decay have a well defined flux, well defined energy 
spectrum, and low backgrounds. Fortunately, accelerators that can provide such a source are 
either being built or are being proposed. The first such accelerator is the 1.4 MW, 1.3 GeV, 
short duty-factor Spallation Neutron Source (SNS) |17j . which is currently being built and will 
be fully commissioned by 2008. The second is a proposal for a 2 MW, 8 GeV, Proton Driver 
(PD) at FNAL Q2|. 

The dominant decay scheme that produces neutrinos from a stopped pion source is 

7r + — ► r = 26 nsec (1) 

followed by 

/i + — > e + i> M f e , r = 2.2 //sec. (2) 

The neutrinos from stopped 7r _ 's are highly suppressed because the negative pions are ab- 
sorbed in the surrounding material. Thus, neutrinos from the ir~ decay chain are significantly 
depleted and can be estimated from the measured z/„, and v c flux. Figure ^ shows the neu- 
trino time and energy spectra from a stopped pion source. As shown in the right hand plot, 
the energy is monoenergetic {E v = 29.8 MeV), while the and u c have known Michel 
decay distributions with an end point energy of 52.8 MeV. Furthermore, the left hand plot 
shows the SNS beam timing and the time distributions for the three neutrino species. One of 
the advantages of the SNS relative to the LAMPF accelerator is the three orders of magnitude 
shorter beam time of 695 nsec (the FNAL PD will have a similar beam time). With a simple 
beam-on timing cut, one can obtain a fairly pure i>u sample, with only a 14% contamination 
of Vu and v e each. This remaining background is easily measured and subtracted. 
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Figure 1: The neutrino time and energy spectra of the different neutrino species produced 
isotropically from a stopped pion source |17j . 
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Table 1: Proton intensities at FNAL and SNS. The numbers are taken from |18) . assuming 
3.15 x 10 7 s/yr operation. 

Table^shows the expected proton rates for both SNS and FNAL beamlines, normalized to 
a full year of running, i.e. 3.15 x 10 seconds. The FNAL Proton Driver proposal is broken down 
into 8 GeV, 2 MW and 5 GeV, 1.25 MW beamlines. The FNAL 8 GeV option will provide about 
2.5 times more protons per year than the SNS. However, the FNAL PD is only a proposal, while 
the SNS is under construction and will be operational by 2008. This makes the SNS a more 
timely option. Furthermore, the SNS is planning for 2014 an upgrade which will deliver 3MW 
to two sources, making the interesting situation of multiple baselines with a single detector. 

A key component of the sterile neutrino measurement is the physical size of the stopped 
pion source, which adds an uncertainty to the neutrino path length. For the SNS, the compact 
liquid mercury target will contribute approximately 25 cm (FWHM) to the neutrino path 
length uncertainty. The FNAL source size should be of similar dimensions to minimize neutrino 
path length uncertainties. 
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4 Active-Sterile Neutrino Oscillation Measurements 



The direct observation of active-sterile neutrino oscillations with Am 2 > 0.1 eV 2 can be 
achieved because of two key features of a short duty-factor stopped pion neutrino source. 
Firstly, all active neutrino species (x = z/ e , v„, v T ) can be efficiently measured using the su- 
perallowed neutral current reaction v x 12 C — ► u x 12 C*(15.11) where the C*(15.11) state 
has a 92 ± 2% branching fraction to the ground state [2U]. A deficit of the neutral current 
rate in 12 C from expectation can only be from oscillations to sterile neutrinos. Secondly, using 
beam timing cuts, we can extract a high purity sample of 29.8 MeV monoenergetic events. 
In (3+1) models with Am 2 > 0.1 eV 2 , active-sterile oscillations can be approximated by two 
neutrino mixing. The oscillation probability can be expressed as: 

-> v„) = 1 - sin 2 26 sin 2 (1.27 Am 2 -J-). (3) 

Given these two features, we can observe active-sterile neutrino oscillations by measuring 
changes in the v x 12 C — » v x 12 C*(15.11) rate as a function of L (in meters). FigureEJshows the 
oscillation length as a function of Am 2 . Thus, a two detector configuration of the appropriate 
size and distance can achieve the desired sensitivity. A near detector greater than three meters 
in length and about 10-20 meters from the source will be sensitive to Am 2 > 20 eV 2 . A second 
large detector (required to increase statistics because of decreasing 1/r 2 neutrino flux) at 
50-100 meters will be sensitive to Am 2 > 0.1 eV 2 . Figure EH shows a simple schematic of a 
detector setup that would allow sensitivity to sterile neutrino oscillations using the v x 12 C — > 
v x 12 C*(15.11) channel. 

A possible detector configuration at the SNS would be first a 25 ton fiducial mass (3.5 m 
on a side) liquid scintillator near detector at 18 m distance from the neutrino source. This is 
the homogeneous detector proposed by the z^-SNS collaboration to measure neutrino carbon 
cross sections |17j . The limited fiducial mass size of the near detector might be improved with 
future developments in flat PMT technology. Next would be a 500 ton fiducial mass liquid 
scintillator (6 m radius) far detector at 60 m distance from the neutrino source. Both designs 
would be similar to an LSND style detector, i.e. a mineral oil detector doped with butyl PBD 
scintillator to increase light output, instrumented with 25% photocathode coverage. This gives 
approximately a 7% energy resolution for reconstruction of the 15.11 MeV gamma ray. Both 
detector masses are quoted in fiducial mass, which is the useful mass available for reliable 
event reconstruction. To estimate NC neutrino rates in the near and far detector, we scale to 
the detection rate calculations in |T2|, which have assumed a z/„ 12 C — > v„ 12 C*(15.11) cross 
section of 2.7 x 10~ 42 cm 2 |19j . The Karmen collaboration has measured this cross section 
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Figure 2: Oscillation length as a function Am 2 where E u = 29.8 MeV. 
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Figure 3: Neutral Current u x 12 C — > u x 12 C* (15.11) disappearance patterns for a two 
detector setup. 
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Table 2: Estimated 12 C -> 12 C*(15.11) events per year at the SNS (1.4 MW) and FNAL (2MW) 
sources, assuming 100% event reconstruction efficiency. 



with limited statistics, getting a value of (3.2 ± 0.5 s t a t ± 0.4 syst ) x 10 -42 cm 2 |21j . consistent 
with the calculated cross section. 

For a proposed experiment at FNAL, we would not have space restrictions, allowing for 
larger near and far detectors. We assume a longer baseline, out to 100 m, and a detector mass 
of 2ktons (8m radius), which would be required to keep roughly the same statistics. This 
would give a fiducial mass of 1.3ktons (7m radius). 

Table |21 shows the number of events expected per year for the u x 12 C — ► u x 12 C*(15.11) 
reaction for the two possible configurations, one at SNS with 1.4 MW, and the other at FNAL 
with 2MW. The event reconstruction efficiency was assumed to be 100%, as the 15.11 MeV 
gamma ray is relatively easy to detect. 

A Monte Carlo calculation was performed to estimate the i/ x 12 C — > u x 12 C*(15.11) oscilla- 
tion sensitivity, which includes smearing of 25 cm (FWHM) from the finite neutrino source size 
and 50 cm (FWHM) from the reconstructed position resolution for the 15.11 MeV gamma ray. 
Also included is background subtraction of charged current (CC) and neutral current (NC) 
events from i?n and v e in the beam window Because we can isolate a clean sample of these 
backgrounds in the time window after the beam, we will have an excellent estimation of their 
number and character, e.g. for the SNS far detector this is about 1300 events per year. This 
background increases the statistical error on the measured v„ 12 C — > 12 C*(15.11) rate by 
at most 45%. The background contribution from cosmic rays is negligible because of the short 
duty factor of the beam. 

Figures 0] and show the (3+1) active-sterile neutrino oscillation sensitivity for a three year 
run at the SNS and FNAL proton driver with two detectors and a 5% normalization systematic 
errors. Both setups achieve a — > u s oscillation sensitivity of 3a for Am 2 > 0.4 eV 2 and 
sin 2 26 > 0.05. This sensitivity is desired if the LSND signal is due to active-sterile neutrino 
oscillations. 

The oscillation sensitivity plots are generated from a simultaneous fit to both the L shape 
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and to an overall event normalization. The L measurement systematic errors are discussed 
above. The event normalization requires knowledge of the flux times neutral current cross 
section (<j> v x onc)- From rate calculations alone, we expect systematic errors on the order of 
10%, which is what LSND was able to achieve jllj . Another technique is to use the charged 
current u c 12 C — ► e~ 12 N gs reaction to normalize the NC rate, as the final state A= 12 nuclear 
form factors are isotopic analogues, leading to systematic errors in the NC/CC ratio of ~ 
3% |21| . Here we assume there are no oscillations of u c . The NC/CC ratio is dominated by 
systematic errors since the event rates for both reactions are large. Besides the theoretical 
ratio error (discussed above), we need to know the relative fiducial volume efficiency and the 
relative detection efficiency of the 15.11 MeV 7 to that of the CC electron plus 12 N beta decay. 
Based on experiences at LSND, these two errors are conservatively estimated at 3% each. 
This gives a total systematic uncertainly in the NC/CC ratio of 5%, which translates into a 
12 C — > z/ju 12 C*(15.11) event rate <j> v x cjnc systematic error of 5%. 

If we consider just a single detector setup at the SNS, either near or far, Figures E] and 
show what might be expected for 5% <f>„ x cjnc errors. Figure El shows that even this one near 
detector has good sensitivity to active-sterile oscillations with parameters of Am 2 > 1 eV 2 
and sin 2 26 > 0.10 at 3a. This interesting measurement can be performed in parallel to the 
neutrino carbon cross section measurement planned for the zz-SNS detector |17j . 

In figures 4 through 9, two black stars are plotted that correspond to the (3+2) model sterile 
neutrino solution Am], = 0.9 eV 2 , sin 2 26 Al = 0.15, Am| x = 22 eV 2 , sin 2 20 5 i = 0.19 0. This 
(3+2) model gives significantly better fits than (3+1) models [I]. Although the sensitivity 
curves are strictly for (3+1) models and cannot be directly applied to the (3+2) solutions, 
it can be seen that all beam/detector combinations will be sensitive to these solutions. It 
should be noted that the favored (3+1) solution from [3j is very close to the lower Am 2 (3+2) 
solution, and so this solution will have coverage at the 3cr to 5a level. A recent (3+1) analysis 
of all the short baseline data, excluding LSND, gives — » v s oscillation plots that exclude the 
upper right regions shown here for Am 2 > leV 2 and sin 2 26 > 0.1 at the 99% level |22j . Hence, 
even ignoring LSND, there is much untested phase space that can be probed for active-sterile 
neutrino oscillations. 

Explicit treatment of (3+2) models is shown in figures Eland El where we present parameter 
determination of the (3+2) solutions. These models cannot be approximated by a two neutrino 
approximation, and consequently sin 2 26 can no longer be interpreted as before. We instead 
define sin 2 26 as 4U 2 n (l — U 2 n ), where U^ n is the relevant element in the neutrino mixing 
matrix and n is the sterile generation label. In these figures, one of the two solutions is 
assumed and fixed. This is a reasonable procedure since the two solutions have quite different 
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characteristics, i.e. the low Am 2 solution is demonstrated by a rate suppression between two 
detectors, while the high Am 2 solution is demonstrated by rapid wiggles in L. This allows 
for easy separation of the two solutions. The resulting plot, therefore, shows the parameter 
determination at a single point in the full four dimensional parameter space, rather than a 
projection in this space. Both figures show that the parameter determination is excellent. 

5 Conclusion 

In the coming year, new physics beyond the Standard Model will be implied if MiniBooNE 
confirms the LSND oscillation anomaly. To understand whether active-sterile neutrino oscil- 
lations are responsible for this new physics will require a new class of oscillation experiments. 
A stopped pion source such as the SNS, coupled with two well positioned large liquid scin- 
tillator detectors, can achieve the desired sensitivity. This is made possible by the fortuitous 
combination of a monoenergetic flux and the efficiently reconstructed neutral current inter- 
action u x 12 C — > z/ x 12 C*(15.11). The SNS two detector setup with 5% <\> v x ctnc systematic 
errors can explore the — > u a oscillation region Am 2 > 0.4 eV 2 and sin 2 26 > 0.05 at the 3a 
level. This would cover the Vn — > u s oscillation solutions indicated by LSND and the (3+2) 
model. Such a measurement would be the first clear indication of the existence of light sterile 
neutrinos, and would herald a new era in neutrino physics. 
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Figure 4: Active-Sterile neutrino oscillation sensitivity with the NC reaction v x 12 C — ► 
u x 12 C*(15.11) for SNS with two detectors, 5% <\> v x onc systematic errors, and three years 
of running. The two black stars correspond to the (3+2) solutions (see text). 




Figure 5: Active-Sterile neutrino oscillation sensitivity with the NC reaction v x 12 C — ► 
u x 12 C*(15.11) for FNAL PD with two detectors, 5% <p u x a^c systematic errors, and 
three years of running. The two black stars correspond to the (3+2) solutions (see text). 
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Figure 6: Active-Sterile neutrino oscillation sensitivity with the NC reaction v x 12 C — > 
u x 12 C*(15.11) for SNS with a near detector, 5% <j> v x cjnc systematic errors, and three 
years of running. The two black stars correspond to the (3+2) solutions (see text). 




Figure 7: Active-Sterile neutrino oscillation sensitivity with the NC reaction v x 12 C — ► 
u x 12 C*(15.11) for SNS with a far detector, 5% <\> v x cjnc systematic errors, and three years 
of running. The two black stars correspond to the (3+2) solutions (see text). 
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Figure 8: The (3+2) active-sterile neutrino oscillation solution contours for SNS with two 
detectors, 5% <p u x a^c systematic errors, and three years of running. See text for definition 
of sin 2 29. 
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Figure 9: The (3+2) active-sterile neutrino oscillation solution contours for FNAL with 
two detectors, 5% tp v x cjnc systematic errors, and three years of running. See text for 
definition of sin 2 26. 
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